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Main  Scientific  Objective 


•Extraction  the  beam-averaged  modal  attenuation  coefficient  5„ 

How  to  do? 

Step  1:  Perform  vertical  plane  wave  beamforming  with  VLA  acoustic 
data 

Step2:  Estimate  the  beamformer  output  spectrum 

Step3:  Calculate  propagation  loss  of  two  different  ranges 

Step4:  Extract  the  beam-averaged  modal  attenuation  coefficient  from 
propagation  loss 

Why  to  do? 

Be  important  for  bottom  back-scattering  matrix  extraction  from 
reverberation  data  as  well  as  propagation  modeling 


Theory  Frame 


The  complex  pressure  in  a  range-independent  waveguide  can  be  written  as: 


A^-l 


p(r,Zj,z„(£>)^Y, iZj )(p„ {z, ) exp[- i{m-kj+Ti ! A)-bj] 


^=0 


For  an  ideal  spatial  filter  response,  the  beamformer  output  can  be  written  as: 


m+Am/  2 


b(r,  Y,  ,  CO)  =  Y,  (r)9„  (Zs )  exp[-  i(cot  -  k„r  +  tc/  4)-  5„r] 


m-Am/2 


The  vertical  beamformer  output  spectrum  is 


I  B(r,  Y„ ,  z, ,  CO) = (r)(p^  (z  J  exp(-25„r) 


n 


a„(r)a^(r)cp„(z^)(p^(z,)exp(-(5„  +6„)r)cos{Ak^r) 


n  m 
n^m 

B1+B2 


of  Oc« 


Theory  Fra^ 


(continued) 


If  B2«Bp  the  beamformer  output  spectrum  is  written  as: 


B(r,  y„ ,  ,  CO)  =  ^  (r)(p2  (z^ )  exp(  -25„r) 


n 


We  define  the  beam-averaged  modal  attenuation  as  follow: 


m+Am/ 2 


'm+Am/2 


exp(-25^r)  =  (r)(^^  (z^ )  exp(-25„r)  /  (r)^^  (z, ) 


m-Am/ 2 


m-Am/ 2 
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Numerical  Simulation 


I 

\ 


of  Oo 


in 


We  consider  a  Pekeris  waveguide  :  H=103m,  CQ=1519m/s, 
Cb=1623m/s,  ?b=1.72,  ab=0.49  dB/m  kHz,  f=1390Hz.  The  103m 

VLA  is  covered  all  depth  of  water  column  with  Im  space  of 
adjacent  element. 


Table!  Theoretical  value  and  calculated  value  of  5^^  in  Pekeris  waveguide  (f=  1390Hz) 


group 

Center 

grazing 

angle 

Group  of 
modes 

Theory  value  of 

s„ 

Calculated  value  of  8^ 

r  1 = 1 0km,r2=20km, 
zs=20m 

r  1 = 1 0km,r2=20km, 
zs=40m 

1 

2.3  i  a 

1-8 

3.367554e-6 

-5.030152e-6 

3.158571e-6 

2 

6.8;  a 

9-16 

2.051041e-5 

3.249568e-5 

1.507090e-5 

3 

11.5;  a 

17-24 

5.064872e-5 

3.641855e-5 

5.898481e-5 

4 

16.2;  a 

25-32 

9.020034e-5 

9.085626e-5 

8.014695e-5 

Due  to  the  effect  of  intermode  interference  ,  8^  of  group  1  can  not  be  estimated  correctly, 
and  the  precisions  of  estimated  5^  of  other  groups  are  poor. 
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Numerical  SimulatioiT 


(continued) 


How  can  we  reduce  the  effects  of  intermode  interference?  One  method  is 
averaged  the  beamformer  output  spectrum  with  f  in  a  narrow  band  to 
smooth  second  item  B2. 


Fig.  3  Effect  of  intermode 
interference  (f= 1 340-  1440Hz, 
z  =20m) 


Fig. 4  Effect  of  intermode 
interference  (f= 1 340-  1440Hz, 
z  =40m) 


Key 


Table  2  Theory  value  and  calculated  value  in  Pekeris  waveguide  (f=  1340- 1440Hz) 


group 

Center 

grazing 

angle 

Group 

of 

modes 

Theory  value 
of8„ 

Calculated  value  oS^^ 

rl=5km, 

r2=10km, 

zs=40m 

rl=10km, 

r2=20km, 

zs=40m 

rl=10km, 

r2=20km, 

zs=20m 

1 

2.3;  a 

1-8 

3.143104e-6 

7.894265e-6 

1.705594e-6 

3.646730e-6 

2 

6.8;  a 

9-16 

1.918959e-5 

1.920825e-5 

2.095787e-5 

1.802004e-5 

3 

11.5i  a 

17-24 

4.754961e-5 

4.799638e-5 

4.731562e-5 

4.815836e-5 

4 

16.2;  a 

25-32 

8.491 166e-5 

8.500221e-5 

8.566655e-5 

8.432893e-5 

Compared  to  Table  1,  the  attenuation  coefficients  of  all  groups  can 
be  estimated,  the  precision  is  also  improved  effectively. 


Numerical  Simulation 


of  0c« 


(continued) 


All  the  calculated  are  larger  than  the 
theoretical  values.  This  is  due  to: 

(1)  effect  of  inter  mode  interference 

(2)  the  short  VLA  covered  the  upper  half 
part  of  the  water  column  is  not  able  to 
capture  the  lower  modes  properly  and 
affected  by  higher  modes  by  leakage. 


Table  3  Theoretical  value  and  calculated  value  of  5^^  in  simulated  shallow  waveguide  (f=1340- 
1440Hz) 


group 

Center 

grazing 

angle 

Group  of 
modes 

Theory  value  of 

8„ 

Calculated  value  of 

r  1 = 1 0km,r2=20km, 
zs=40m,  103mVLA 

r  1 = 1 0km,r2=20km, 
zs=40m,57mVLA 

1 

2.3;  a 

1-8 

6.011310e-6 

2.193981e-5 

1.630264e-5 

2 

6.8;  a 

9-16 

2.272945e-5 

3.075 102e-5 

4.356382e-5 

3 

11.5;  a 

17-24 

3.601585e-5 

5.283101e-5 

6.074829e-5 

4 

16.2;  a 

25-32 

8.115852e-5 

8.317109e-5 

1.071509e-4 
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ECS  Real  data  analysis 


Experiment  setup 


)  Not  used  in  this  paper 


sound  speed(in/s) 
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ECS  Real  data  analysis; 


r 

1  (continued) 


Experiment  setup 

•Receiving  array:  60-element  VLA,  from  3m  to  60m 

•Source:  locating  at  different  depths  Zg=5m,  20m,  40m,  and 
different  ranges  r=5km,  10km,  20km  and  30~40km 

•  Water  depth:  H=103m 

•  Signal  forms:  1 100-2000Hz  PRN,  630Hz  CW,  630Hz  PCW, 
1250-2500HZ  LFM,  and  explosive  sources 

In  our  talk,  only  signals  of  1 100-2000Hz  PRN  at  ranges  of 
5km  and  10km,  depths  of  20m  and  40m  are  used  to  estimated 
beam-averaged  5„  now. 
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ECS  Real  data  analysis 


Fig.6  Vertical  beam  pattern  with  20m  source  depth 


Table  4  The  estimated  value  of  5^^  in  real  shallow  waveguide  (f=1340- 
1440Hz) 


group 

Center 

grazing 

angle 

Group  of 
modes 

Estimated  value  of6j^ 

rl=5km,r2=10km, 

zs=20m 

rl=5km,r2=10km, 

zs=40m 

1 

2.3i  a 

1-8 

1.195583e-4 

2.785835e-5 

2 

6.8;  a 

9-16 

2.447215e-4 

1.224623e-4 

3 

11.5i  a 

17-24 

1.853596e-4 

1.175435e-4 

4 

16.2i  a 

25-32 

2.953209e-4 

1.934406e-4 
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Summary 


•If  the  intermode  sum  can  be  ignored,  the  beamformer  output 
spectrum  of  two  different  ranges  can  be  used  to  extract  beam- 
averaged  modal  attenuation  coefficient. 

•When  vertical  beamformer  is  performed  in  a  narrow  band  not  in 
single  tone,  the  intermode  sum  can  be  reduced  effectively,  as 
proved  by  the  numerical  examples,  but  it  is  still  large,  should  be 
decreased  further. 


•From  numerical  simulation  and  ECS  real  data  analysis  results,  we 
should  extract  beam-averaged  model  attenuation  coefficient  from 
data  collected  by  longer  VLA  covered  all  water  column  in  a 
broader  band. 
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Thank  you 


